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Abstract 

In an inhomogeneous magnetised plasma the transport of energy and particles perpendicular to 
the magnetic field is in general mainly caused by quasi two-dimensional turbulent fluid mixing. The 
physics of turbulence and structure formation is of ubiquitous importance to every magnetically 
confined laboratory plasma for experimental or industrial application. Specifically, high temper- 
ature plasmas for fusion energy research are also dominated by the properties of this turbulent 
transport. Self-organisation of turbulent vortices to mesoscopic structures like zonal flows is re- 
lated to the formation of transport barriers that can significantly enhance the confinement of a 
fusion plasma. This subject of great importance in research is rarely touched on in introductory 
plasma physics or continuum dynamics courses. Here a brief tutorial on 2D fluid and plasma tur- 
bulence is presented as an introduction to the field, appropriate for inclusion in undergraduate and 
graduate courses. 
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I. INTRODUCTION 



Turbulence is a state of spatio-temporal chaotic flow generically attainable for fluids with 
access to a sufficient source of free energy. A result of turbulence is enhanced mixing of the 
fluid which is directed towards a reduction of the free energy. Mixing typically occurs by 
formation of vortex structures on a large range of spatial and temporal scales, that span 
between system, energy injection and dissipation scales [l-3]. 

Fluids comprise the states of matter of liquids, gases and plasmas J4]. A common free 
energy source that can drive turbulence in neutral (or more precisely: non-conducting) 
fluids is a strong enough gradient (or "shear") in flow velocity, which can lead to vortex 
formation by Kelvin-Helmholtz instability. Examples for turbulence occuring from this type 
of instability are forced pipe flows, where a velocity shear layer is developing at the wall 
boundary, or a fast jet streaming into a stationary fluid. Another source of free energy is 
a thermal gradient in connection with an aligned restoring force (as in liquids heated from 
below in a gravity field) that leads to Rayleigh-Benard convection 

Several routes for the transition from laminar flow to turbulence in fluids have been 
proposed. For example, in some specific cases the Ruelle-Takens scenario occurs, where 
by linear instability through a series of a few period doubling bifurcations a final nonlinear 
transition to flow chaos is observed when a control parameter (like the gradient of velocity or 
temperature) is increased |(|. For other scenarios, like in pipe flow, a sudden direct transition 
by subcritical instability to fully developed turbulence or an intermittent transition are 
possible 0,3- 

The complexity of the flow dynamics is considerably enhanced in a plasma compared 
to a non-conducting fluid. A plasma is a macroscopically neutral gas composed of many 
electrically charged particles that is essentially determined by collective degrees of freedom 
0, Q. Space and laboratory plasmas are usually composed of positively charged ions 
and negatively charged electrons that are dynamically coupled by electromagnetic forces. 
Thermodynamic properties are governed by collisional equilibration and conservation laws 
like in non-conducting fluids. The additional long-range collective interaction by spatially 
and temporally varying electric and magnetic fields allows for rich dynamical behaviour of 
plasmas with the possibility for complex flows and structure formation in the presence of 
various additional free energy sources 



2 



FIG. 1: Computation of decaying two-dimensional fluid turbulence, showing contours 
of vorticity u = V x u. 



The basic physics of plasmas in space, laboratory and fusion experiments is introduced 
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y simple, 



in detail in a variety of textbooks (e.g. in Refs. 

Although the dynamical equations for fluid and plasma flows can be conceptual 
they are highly nonlinear and involve an infinite number of degrees of freedom 
alytical solutions are therefore in general impossible. The description of fluid and plasma 
dynamics mainly relies on statistical and numerical methods. 
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il. CONTINUUM DYNAMICAL THEORY OF FLUIDS AND PLASMAS 

Computational models for fluid and plasma dynamics may be broadly classified into three 
major categories: 

• (1) Microscopic models: many body dynamical description by ordinary differential 
equations and laws of motion; 

• (2) Mesoscopic models: statistical physics description (usually by integro-differential 
equations) based on probability theory and stochastic processes; 

• (3) Macroscopic models: continuum description by partial differential equations based 
on conservation laws for the distribution function or its fluid moments. 



Examples of microscopic models are Molecular Dynamics (MD) methods for neutral flu- 
ids that model motion of many particles connected by short range interactio ns 116 1 , or 
Particle- In-Cell (PIC) methods for plasmas including electromagnetic forces 17], [l8|. Such 
methods become important when relevant microscopic effects are not covered by the averag- 
ing procedures used to obtain meso- or macroscopic models, but they usually are intrinsical 
computationally expensive. 

Complications result from multi-particle or multi-scale interactions. Mesoscopic mod- 
elling treats such effects on the dynamical evolution of particles (or modes) by statistical 



assumptions on the interactions 



19]. 



These may be implemented either on the macroscale as spectral fluid closure schemes 
like, for example, in the Direct Interaction Approximation (DIA), or on the microscale as 
advanced collision operators like in Fokker-Planck models. An example of a mesoscopic 
computational model for fluid flow is the Lattice Boltzmann Method (LBM) that combines 
free streaming particle motion by a minimalistic discretisation in velocity space with suitable 
models for collision operators in such a way that fluid motion is recovered on macroscopic 
scales. 

Macroscopic models are based on the continuum description of the kinetic distribution 
function of particles in a fluid, or of its hydrodynamic moments. The continuum modelling 
of fluids and plasmas is introduced in more detail below. 

Computational methods for turbulence simulations have been developed within the frame- 
work of all particle, mesoscopic or continuum models. Each of the models has both advan- 
tages and disadvantages in their practical numerical application. The continuum approach 
can be used in situations where discrete particle effects on turbulent convection processes are 
negligible. This is to some approximation also the case for many situations and regimes of 
interest in fusion plasma experiments that are dominated by turbulent convective transport, 
in particular at the (more collisional) plasma edge. 

Within the field of Computational Fluid Dynamics the longest experience and broadest 
applications have been obtained with continuum methods {20 1. Many numerical continuum 
methods that were originally developed for neutral fluid simulation have been straightfor- 
wardly applied to plasma physics problems 21]. 

In continuum kinetics, the time evolution of the single-particle probability distribution 
function /(x, v,t) for particles of each species (e.g. electrons and ions in a plasma) in 
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the presence of a mean force field F(x, t) and within the binary col 
(modelled by an operator C) is described by the Boltzmann equation 



ision approximation 
22 1 



{d t + vd x + F-d v )f = Cf. (1) 

In a plasma the force field has to be self-consistently determined by solution of the Maxwell 
equations. Usually, kinetic theory and computation for gas and plasma dynamics make 
use of further simplifying approximations that considerably reduce the complexity: in the 
Vlasov equation binary collisions are neglected (C = 0), and in the drift-kinetic or gyro- 
kinetic plasma equations further reducing assumptions are taken about the time and space 
scales under consideration. 

The continuum description is further simplified when the fluid can be assumed to be in 
local thermodynamic equilibrium. Then a hierarchical set of hydrodynamic conservation 



equations is obtained by construction of moments over velocity space 23|, |24j. In lowest 
orders of the infinite hierarchy, the conservation equations for mass density n(x, t), momen- 
tum nu(x, t) and energy density £(x, t) are obtained. Any truncation of the hierarchy of 
moments requires the use of a closure scheme that relates quantities depending on higher 
order moments by a constitutive relation to the lower order field quantities. 

An example of a continuum model for neutral fluid flow are the Navier-Stokes equations. 
In their most widely used form (in particular for technical and engineering applications) the 
assumptions of incompressible divergence free flow (i.e., n is constant on particle paths) and 



of an isothermal equation of state are taken 25]. 



Then the description of fluid flow can be reduced to the solution of the (momentum) 
Navier-Stokes equation 

(d t + u • V) u = - VP + z/Au (2) 

under the constraints given by 

V ■ u = (3) 

and by boundary conditions. Most numerical schemes for the Navier-Stokes equation require 
solution of a Poisson type equation for the normalised scalar pressure P = p/po in order to 
guarantee divergence free flow. 

The character of solutions for the Navier-Stokes equation intrinsically depends on the 
ratio between the dissipation time scale (determined by the kinematic viscosity v) and the 
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mean flow time scale (determined by the system size L and mean velocity U), specified by 
the Reynolds number 

Re = LU/u. (4) 

For small values of R e the viscosity will dominate the time evolution of u(x) in the Navier- 
Stokes equation, and the flow is laminar. For higher R e the advective nonlinearity is dom- 
inant and the flow can become turbulent. The Rayleigh number has a similar role for the 
onset of thermal convective turbulence I™. 



III. DRIFT-REDUCED TWO-FLUID EQUATIONS FOR PLASMA DYNAMICS 

Flow instabilities as a cause for turbulence, like those driven by flow shear or thermal con- 
vection, do in principle also exist in plasmas similar to neutral fluids |26|], but are in general 
found to be less dominant in strongly magnetised plasmas. The most important mechanism 
which results in turbulent transport and enhanced mixing relevant to confinement in mag- 



netised plasmas [27H29l| is an unstable growth of coupled wave-like perturbations in plasma 
pressure p and electric fields E. The electric field forces a flow with the ExB ( "E-cross-B" ) 
drift velocity 



v E xB = -^E x B (5) 
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& 

of the plasma perpendicular to the magnetic field B. A phase shift, caused by any inhibition 
of a fast parallel Boltzmann response of the electrons, between pressure and electric field 
perturbation in the presence of a pressure gradient can lead to an effective transport of 
plasma across the magnetic field and to an unstable growth of the perturbation amplitude. 
Nonlinear self-advection of the ExB flow and coupling between perturbation modes ("drift 
waves" ) can finally lead to a fully developed turbulent state with strongly enhanced mixing. 

A generic source of free energy for magnetised laboratory plasma turbulence resides in 
the pressure gradient: in the core of a magnetic confinement region both plasma density 
and temperature are usually much larger than near the bounding material wall, resulting 
in a pressure gradient directed inwards to the plasma center. Instabilities which tap this 
free energy tend to lead to enhanced mixing and transport of energy and particles down 



the gradient 30]. For magnetically confined fusion plasmas, this turbulent convection by 



ExB drift waves often dominates collisional diffusive transport mechanisms by orders of 



magnitude, and essentially determines energy and particle confinement properties 
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32] 



The drift wave turbulence is regulated by formation of mesoscopic streamers and zonal 
structures out of the turbulent flows 



Continuum models for drift wave turbulence have to capture the different dynamics of 
electrons and ions parallel and perpendicular to the magnetic field and the coupling between 
both species by electric and magnetic interactions [34, 35]. Therefore, a single-fluid magneto- 



hydrodynamic (MHD) model can not appropriately describe drift wave dynamics: one has 
to refer to a set of two-fluid equations, treating electrons and ions as separate species, 
although the plasma on macroscopic scales remains quasi-neutral with nearly identical ion 
and electron density, rii n e = n. 

The two-fluid equations require quantities like collisional momentum exchange rate, pres- 
sure tensor and heat flux to be expressed by hydrodynamic moments based on solution of a 
kinetic (Fokker- Planck) model. The most widely used set of such fluid equations has been 



derived by Braginskii [36] and is e.g. presented in brief in Ref. [37 ]. 

The most general continuum descriptions for the plasma species, based either on the ki- 
netic Boltzmann equation or on the hydrodynamic moment approach like in the Braginskii 
equations, are covering all time and space scales, including detailed gyro-motion of particles 



around the magnetic field lines, and the fast plasma 
mental observation it is on the other hand evident 38 



requency oscillations. From experi- 



42j, that the dominant contributions 



to turbulence and transport in magnetised plasmas originate from time and space scales 
that are associated with frequencies in the order the drift frequency u ~ (p s /L±)Qi, that 
are much lower than the ion gyro-frequency Qi = qiB/Mi by the ratio between drift scale 



p s = y/T e Mi/ (eB) to gradient length L±: 

u ~ d t < Vti < Vt e . (6) 

Under these assumptions one can apply a "drift ordering" based on the smallness of the 
order parameter 5 = < 1. This can be introduced either on the kinetic level, resulting 
in the drift-kinetic model, or on the level of two-fluid moment equations for the plasma, 



resu 
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ting in the "drift-reduced two-fluid equations" , or simply called "drift wave equations" 
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44j: neglect of terms scaling with 5 in higher powers than 2 considerably simplifies 
both the numerical and analytical treatment of the dynamics, while retaining all effects of 
the perpendicular drift motion of guiding centers and nonlinear couplings that are necessary 
to describe drift wave turbulence. 
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For finite ion temperature, the ion gyro-radius pi = ^TiMi/(eB) can be of the same 
magnitude as typical fluctuation scales, with wave numbers found around k±p s ~ 0.3 in the 
order of the drift scale 

p » = -Jb~- (7) 

Although the gyro-motion is still fast compared to turbulent time scales, the ion orbit then 
is of similar size as spatial variations of the fluctuating electrostatic potential. Finite gyro- 
radius (or "finite Larmor radius", FLR) effects are captured by appropriate averaging pro- 
cedures over the gyrating particle trajectory and modification of the polarisation equation, 
resulting in "gyrokinetic" or "gyrofluid models" for the plasma. 

IV. TURBULENT VORTICES AND MEAN FLOWS 

The prevalent picture of drift wave turbulence is that of small-scale, low-frequency ExB 
vortices in the size of several gyro-radii, that determine mixing and transport of the plasma 
perpendicular to the magnetic field across these scales. 

Beyond that, turbulence in magnetised plasmas exhibits large-scale structure formation 
that is linked to this small-scale eddy motion: The genesis of mean zonal flow structures 
out of an initially nearly homogeneous isotropic vortex field and the resulting shear-flow 
suppression of the driving turbulence is a particular example of a self-organising regulation 



process in a dynamical system |27l . 145m lj . The scale of these macroscopic turbulent zonal 
flows is that of the system size, setting up a radially localised differential ExB rotation of 
the whole plasma on an entire toroidal flux surface. 

Moreover, the process of self-organisation to zonal flow structures is thought to be a 
main ingredient in the still unresolved nature of the L-H transition in magnetically confined 



toroidal plasmas for fusion research 52]. The L-H transition is experimentally found to be 
a sudden improvement in the global energy content of the fusion plasma from a low to high 
(L-H) confinement state when the central heat source power is increased above a certain 



threshold 
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-561]. The prospect of operation in a high confinement H-mode is one of the 
main requirements for successful performance of fusion experiments like ITER. 

The mechanism for spin-up of zonal flows in drift wave turbulence is a result of the 
quasi two-dimensional nature of the nonlinear ExB dynamics, in connection with the double 
periodicity and parallel coupling in a toroidal magnetised plasma. 
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FIG. 2: Generation of mean sheared Bows from drift wave turbulence within the poloidal cross- 
section of a magnetised plasma torus. 

Basic concepts and terminology for the interaction between vortex turbulence and mean 
flows have been first developed in the context of neutral fluids. It is therefore instructive to 
briefly review these relations in the framework of the Navier-Stokes Eq. fl2]) before applying 
them to plasma dynamics. 

Small (space and/or time) scale vortices and mean flows may be separated formally by 
an ansatz known as Reynolds decomposition, 



u = U + u, 



(8) 



splitting the flow velocity into a mean part U = (u), averaged over the separation scale, 
and small-scale fluctuations u with (u) = 0. While the averaging procedure, (...), is math- 
ematically most unambiguous for the ensemble average, the physical interpretation in fluid 
dynamics makes a time or space decomposition more appropriate. Applying this averaging 
on the Navier-Stokes Eq. ([2]), one obtains the Reynolds equation (or: Reynolds averaged 
Navier-Stokes equation, RANS): 



(dt + U • V) U = - VP + VR + z/AU 



(9) 



This mean flow equation has the same structure as the original Navier-Stokes equation with 
one additional term including the Reynolds stress tensor Rij = (iiiUj). Momentum transport 



between turbulence and mean flows can thus be caused by a mean pressure gradient, viscous 
forces, and Reynolds stress. A practical application of the RANS is in Large Eddy Simulation 
(LES) of fluid turbulence, which efficiently reduces the time and space scales necessary for 
computation by modelling the Reynolds stress tensor for the smaller scales as a local function 
of the large scale flow. LES is however not applicable for drift wave turbulence computations, 
as here in any case all scales down to the effective gyro-radius (or drift scale p s ) have to be 
resolved in Direct Numerical Simulation (DNS). 



V. TWO-DIMENSIONAL FLUID TURBULENCE 



Turbulent flows are generically three-dimensional. In some particular situations the de- 
pendence of the convective flow dynamics on one of the Cartesian directions can be negligible 
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Exam- 



compared to the others and the turbulence becomes quasi two-dimensional 
pies for such 2D fluid systems are thin films (e.g. soap films), rapidly rotating stratified 
fluids, or geophysical flows of ocean currents and the (thin) planetary atmosphere. In par- 
ticular, also the perpendicular ExB dynamics in magnetised plasmas behaves similar to a 
2D fluid {59]. The two-dimensional approximation of fluid dynamics not only simplifies the 
treatment, but moreover introduces distinctly different behaviour. 
The major difference can be discerned by introducing the vorticity 

u> = V x u (10) 

and taking the curl of the Navier-Stokes Eq. ([2]) to get the vorticity equation 

{dt + u • V) u = (w • V) u + uAuj. (11) 

In a two-dimensional fluid with v = v x e x + v y e y and v z — the vorticity reduces to u) = we z 
with w = d x v y — d y v x . The vortex stretching and twisting term (u> • V) u is zero in 2D, thus 
eliminating a characteristic feature of 3D turbulence. For unforced inviscid 2D fluids then 
due to (dt + u • V) uj = the vorticity w is constant in flows along the fluid element. This 
implies conservation of total enstrophy W = J(l/2)|oj| 2 (ix in addition to the conservation 
of kinetic flow energy E = J (l/2)|u| 2 <ix. 

The 2D vorticity equation can be further rewritten in terms of a scalar stream function 
(f) that is defined by (v x ,v y ) = (d y , —d x )<ft so that w = V 2 0, to obtain 

d t w + [cp, w] = uAw. (12) 
10 



Here the Poisson bracket [a, b] = d y a d x b — d x a d y b is introduced. For force driven flows a 
term given by the curl of the force adds to the right hand side of Eq. f fT2|) . Although the 
pressure is effectively eliminated from that equation, it is still necessary to similarly solve a 
(nonlocal) Poisson equation for the stream function. For ExB flows in magnetised plasmas 
the stream function is actually represented by the fluctuating electrostatic potential. 

The energetics of homogeneous 3D turbulence is usually understood in terms of a direct 
cascade of energy from the injection scale down to small (molecular) dissipation scales 
large vortices break up into smaller ones due to mutual stretching and shearing. In terms of 
the Reynolds Eq. ((9]) this means that the Reynolds stress transfer is usually negative, taking 
energy out of mean flows into small scale vortices. The interaction between scales takes 
place basically by three-mode coupling maintained by the convective quadratic nonlinearity. 
This leads to the generic Kolmogorov k~ 5 ^ 3 power spectrum of Fourier components E{k) = 
J dx (l/2)|u(x)| 2 exp(— ikx) in the cascade range of 3D turbulence when energy injection 
and dissipation scales are well separated by a high Reynolds number 60]. 

In two dimensions the behaviour is somewhat different: Kraichnan, Leith and Batchelor 
61 6s| have conjectured that the energy has an inverse cascade property to scales larger 



than the injection. Smaller vortex structures self-organise to merge into bigger ones as 
a result of the absence of vortex stretching. For unforced turbulence the Reynolds stress 
transfer is on the average positive and into the mean flows. 

The classical theory of 2D fluid turbulence by Kraichnan et al. predicts a k~ 3 energy 
spectrum and a k~ x enstrophy spectrum in the inertial range. Numerical simulations of 2D 
Navier-Stokes turbulence however rather find, for example, a fc _5//3 inverse cascade for energy 
on large scales and a k~ 3 direct cascade for enstrophy on small scales jfi^j], modifying the 
classical predictions due to the existence of intermittency and coherent structures, although 
the extent of the modification is still under discussion. The (limiting or periodic) domain 
boundary in 2D simulations has also been found to have stronger influence than for the 3D 
case. 

Periodicity in one dimension of the 2D system can lead to the spin up of sustained zonal 
structures of the mean flow out of the turbulence by inverse cascade. Prominent examples 
of zonal flows in planetary atmospheric dynamics are the well visible structures spanning 
around the planet Jupiter approximately along constant latitude, and jet streams in the 
earth's atmosphere. Zonal flows are also observed in fluids rotating in a circular basin. 
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FIG. 3: Example for a fluid simulation of 2D grid turbulence with a Lattice-Boltzmann code \6a] : 
a high Reynolds number Bow with R e = 5000 is entering from the left of the domain and passes 
around a grid of obstacles. The shading depicts vorticity. In the near field directly behind the 
grid the particular vortex streets can be distinguished. In the middle of the domain neighbouring 
eddies are strongly coupled to a quasi homogeneous (statistically in the perpendicular direction ) 
vortex held. On the far right side eddies decay into larger structures in the characteristic way of 
2D turbulence. The simulation agrees well with Bowing soap film experiments [6a] . 



Drift wave turbulence in magnetised plasmas also has basically a 2D character and ex- 
hibits zonal structure formation in the poloidally and toroidally periodic domain on magnetic 
flux surfaces of a torus. These zonal plasma flows have finite radial extension and constitute 
a differential, sheared rotation of the whole plasma on flux surfaces. 



VI. TURBULENCE IN MAGNETISED PLASMAS 

Drift wave turbulence is nonlinear, non-periodic motion involving disturbances on a back- 
ground thermal gradient of a magnetised plasma and eddies of fluid like motion in which 
the advecting velocity of all charged species is the ExB velocity. The disturbances in the 
electric field E implied by the presence of these eddies are caused by the tendency of the 
electron dynamics to establish a force balance along the magnetic field B. 

Pressure disturbances have their parallel gradients balanced by a parallel electric field, 
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whose static part is given by the parallel gradient of the electrostatic potential. This po- 
tential in turn is the stream function for the ExB velocity in drift planes, which are locally 
perpendicular to the magnetic field. The turbulence is driven by the background gradient, 
and the electron pressure and electrostatic potential are coupled together through parallel 
currents. Departures from the static force balance are mediated primarily through elec- 
tromagnetic induction and resistive friction, but also the electrons inertia, which is not 
negligible. 

The dynamical character of cross-field ExB drift wave turbulence in the edge region of 
a tokamak plasma is governed by this electromagnetic and dissipative effects in the parallel 
response. 

The most basic drift-Alfven (DALF) model to capture the drift wave dynamics includes 
nonlinear evolution equations of three fluctuating fields: the electrostatic potential <ft, elec- 
tromagnetic potential A\\, and density n. The tokamak edge usually features a more or less 
pronounced density pedestal, and the dominant contribution to the free energy drive to the 
turbulence by the inhomogeneous pressure background is thus due to the density gradient. 

On the other hand, a steep enough ion temperature gradient (ITG) does not only change 
the turbulent transport quantitatively, but adds new interchange physics into the dynamics. 
In addition, more field quantities have to be treated: parallel and perpendicular temperatures 
Tii and Tj_ and the associated parallel heat fluxes, for a total of six moment variables for each 
species. Finite Larmor radius effects introduced by warm ions require a gyrofluid description 
of the turbulence equations. 

Both the resistive DALF and the ITG models can be covered by using the six-moment 
electromagnetic gyrofluid model GEM by Scott [671 ]. but for basic studies it is also widely 
used in its more economical two-moment version for scenarios where the DALF model is ap- 
plicable {(38 1 . The gyrofluid model is based upon a moment approximation of the underlying 
gyrokinetic equation. 



The first comp 
Dorland et al. 



ete six- moment gyrofluid formulation was given for slab geometry 



and later extended by Beer et al. to incorporate toroidal effects 70] 



by 



using a ballooning-based form of flux surface geometry 



7l|. 



Electromagnetic induction and electron collisionality were then included to form a more 



general gyrofluid for edge turbulence by Scott 



721 ]. with the geometry correspondingly re- 



placed by the version from the edge turbulence work, which does not make ballooning 
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core edge SOL 




FIG. 4: Left: A cut of a torus shows the poloidal cross section with minor radius r, major radius R, 
poloidal angle 6 and toroidal angle £. Right: typical radial profiles of density n(r) and temperature 
T(r) in the toroidal plasma of a tokamak. The density often is nearly constant in the "core" region 
and shows a pronounced steep-gradient pedestal in the plasma "edge" region. The outer scrape-off 
layer ("SOL") connects the plasma with materials walls. 



assumptions and in particular represents slab and toroidal mode types equally well and does 
not require radial periodicit y |73l| . Energy conservation considerations were solidified first 

and recently for the six-moment version in Ref. 



for the two-moment version 168 



67] 



VII. BASIC DRIFT WAVE INSTABILITY 



Destabilization of the ExB drift waves occurs when the parallel electron dynamics deviates 
from a fast "adiabatic" response to potential perturbations, resulting in a phase shift between 
the density and potential fluctuations. In this section the basic linear instability mechanism 
is discussed in the most basic electrostatic, cold ion limit for a straight magnetic field. 

Figure schematically shows a localized perturbation of plasma pressure p (left) that 
results in a positive potential perturbation > (middle) due to ambipolar diffusion. For 
typical tokamak parameters it is found that the perturbation scale A ^> = \/eqTQ (ne 2 ) 
is much larger than the Debye length X D , so that quasi neutrality n« ~ n e = n can be 
assumed. 

In accordance with the stationary parallel electron momentum balance equation 

-en E N - V|,p e = (13) 
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with p e = h e T e , the isothermal electrons try to locally establish along the field line a Boltz- 
mann relation n e = no(r) exp(e0/T e ). Under quasi neutrality ni=n e =rio(r)+n e , where uq(t) 
is the (in general radially varying) background density. 

Without restrictions on the parallel electron dynamics (like e.g. due to collisions, Alfven 
waves or kinetic effects like Landau damping and particles trapped in magnetic field inho- 
mogeneities) this balance is established instantaneously on the drift time scale and is usually 
termed an "adiabatic response". 

Already at homogeneous background density the perturbation convects the plasma with 
the ExB drift velocity vi = wexB — (-B~ 2 )E x B equal for electrons and ions. When a per- 
pendicular background pressure gradient is present, the perturbed structure propagates 
in the electron diamagnetic drift direction ~ x B. 

In the continuity equation 

d t n + V • (nv) = (14) 

for cold ions in a homogeneous magnetic field and by neglecting ion inertia the only contri- 
bution to the velocity is the perpendicular ExB drift velocity v# = — B~ 2 (V<f) x B). Using 
the Boltzmann relation in Eq. (1141) one gets 



d t n exp I — | - V • 



x B) n exp ( ^ 



0. 



(15) 



and due to the straight B = Be» it is obtained 

T, 



eB 



(d r lnn Q ) d e <p = 0. 



(16) 



Assuming a perturbation periodical in the electron diamagnetic drift coordinate 9 with 
= $ exp[— iut + iked], the electron drift wave frequency is found to be 



w* e = ^-^-h = y-[p s k 6 ] = Y-£li[p a k$\. 



(17) 



Here the density gradient length L n = (d r ln n Q )~ l and the drift scale p s = y/rriiT e / (eB), 
representing an ion radius at electron temperature, have been introduced. 

The motion of the perturbed structure perpendicular to magnetic field and pressure 
gradient in the electron diamagnetic drift direction kg ~ Vp x B is in this approximation 
still stable and does so far not cause any transport down the gradient. 
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FIG. 5: Basic drift waves mechanism: (1) an initial pressure perturbation p leads to an ambipolar 
loss of electrons along the magnetic field B, whereas ions remain more immobile. (2) The resulting 
electric field E = — V</> convects the whole plasma with vexB around the perturbation in the plane 
perpendicular to B. (3) In the presence of a pressure gradient, p propagates in electron diamagnetic 
drift direction with v* ~ Vp x B. This "drift wave" is stable if the electrons establish <fi according 
to the Boltzmann relation without delay ("adiabatic response"). A non-adiabatic response due 
to collisions, magnetic flutter or wave-kinetic effects causes a phase shift between p and <j>. The 
ExB velocity is then effectively transporting plasma down the gradient, enhances the principal 
perturbation and leads to an unstable growth of the wave amplitude 

The drift wave is destabilized only when a phase shift 5^ between potential and density 
perturbation is introduced by "non-adiabatic" electron dynamics 

e<f) 

n e = n (l - idk) =- (18) 

1 e 

The imaginary term iS^ in general is an anti-hermitian operator and describes dissipation of 
the electrons, that causes the density perturbations to proceed the potential perturbations 
in 6 by slowing down the parallel equilibration. This leads to an exponential growth of the 
perturbation amplitude by exp(7 fc t) with linear growth rate 7 fc ~ (Wfc- 

Parallel electron motion also couples drift waves to shear Alfven waves, which are parallel 
propagating perpendicular magnetic field perturbations. With the vector potential Au as a 
further dynamic variable, the parallel electric field E\\, parallel electron motion, and nonlin- 
early the parallel gradient are modified. The resulting nonlinear drift- Alfven equations are 
discussed in the following section. 

The stability and characteristics of drift waves and resulting plasma turbulence are further 
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influenced by inhomogeneities in the magnetic field, in particular by field line curvature and 
shear. The normal and geodesic components of field line curvature have different roles 
for drift wave turbulence instabilities and saturation. The field gradient force associated 
with the normal curvature, if aligned with the plasma pressure gradient, can either act to 
restore or amplify pressure gradient driven instabilities by compression of the fluid drifts, 
depending on the sign of alignment. The geodesic curvature describes the compression of 
the field strength in perpendicular direction on a flux surface and is consequently related to 
the compression of large-scale (zonal) ExB flows. 

Transition from stable drift waves to turbulence has been studied experimentally in linear 
and simple toroidal magnetic field configurations, and by direct numerical simulation. 

Experimental investigations in a magnetized low-beta plasma with clindrical geometry 
by Klinger et al. have demonstrated that the spatiotemporal dynamics of interacting desta- 
bilised travelling drift wave follows a bifurcation sequence towards weakly developed tur- 
bulence according to the Ruelle-Takens-Newhouse scenario jz^. The relationship between 
observations made in linear magnetic geometry, purely toroidal geometry and magnetic con- 
finement is discussed in Ref. 75j, where the role of large-scale fluctuation structures has 
been highlighted. The role of parallel electron dynamics and Alfven waves for coherent drift 
modes and drift wave turbulence have been studied in a collisionality dominated high-density 
helicon plasma [76]. Measurements of the phase coupling between spectral components of 
interchange unstable drift waves at different frequencies in a basic toroidal magnetic field 
configuration have indicated that the transition from a coherent to a turbulent spectrum is 

n 

mainly due to three- wave interaction processes [77j. 

The competition between drift wave and interchange physics in ExB drift turbulence has 
been studied computational^ in tokamak geometry with respect to the linear and nonlinear 
mode structure by Scott [78]. A quite remarkable aspect of fully developed drift wave tur- 
bulence in a sheared magnetic field lying in closed surfaces is its strong nonlinear character, 
which can be self-sustaining even in the absence of linear instabilities 79] . This situation of 
self- sustained plasma turbulence does not have any analogy in neutral fluid dynamics and, as 
shown in numerical simulations by Scott, is mostly applicable to tokamak edge turbulence, 
where linear forcing is low enough so that the nonlinear physics can efficiently operate 
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VIII. DRIFT- ALFVEN TURBULENCE SIMULATIONS FOR FUSION PLASMAS 



The model DALF3 by Scott [80|], in the cold ion approximation without gyrofluid FLR 
corrections, represents the four field version of the dissipative drift-Alfven equations, with 
disturbances (denoted by the tilde) in the ExB vorticity Cl, electron pressure p e , parallel 
current Jn, and parallel ion velocity u» as dependent variables. The equations are derived 
under gyro/drift ordering, in a three dimensional globally consistent flux tube geometry 



m 



8l|, and appear (in cgs units as used in the references) as 
nMr 2 

{dt +VE . V)n = V||J|| _ /Cfe); (19) 

| _ fjl _ j ~ ~ 

-d t A\\ + ^(d t + v E -V)]\\ = — V||(p e + p e )-V||0-77||i||, (20) 

(d t + V E ■ V) (p e + p e ) = "^V|| j\\ - PeV||U|| + Pe^{4>) ~ "^^(Pe), 

(21) 

nM t (d t + v E ■ V) U|, = -V||(p e +p e ), (22) 

with the parallel magnetic potential An given by ju = —(c/Att)'V 2 ± A\\ through Ampere's 
law, and the vorticity Q = V 2 ± 4>. Here, 7711 is the Braginskii parallel resistivity, m e and Mi 
are the electron and ion masses, n is the electron (and ion) density, and T e is the electron 
temperature with pressure p e = nT e . The dynamical character of the system is further 
determined by a set of parameters characterising the relative role of dissipative, inertial and 
electromagnetic effects in addition to the driving by gradients of density and temperature. 

The flux surface geometry of a tokamak enters into the fluid and gyrofluid equations 
via the curvilinear generalisation of differentiation operators and via inhomogeneity of the 
magnetic field strength B. The different scales of equilibrium and fluctuations parallel and 
perpendicular to the magnetic field motivate the use of field aligned flux coordinates. The 
differential operators in the field aligned frame are the parallel gradient 

V|| = (1/£)(B + Bl)-V, (23) 

with magnetic field disturbances Bj_ = (— 1/S)B x VAn as additional nonlinearities, the 
perpendicular Laplacian 

Vl = V-[(-l/5 2 )Bx(BxV)], (24) 
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and the curvature operator 

K, = V- {{c/B 2 )B x V)]. (25) 

The DALF equations constitute the most basic model containing the principal interac- 
tions of dissipative drift wave physics in a general closed magnetic flux surface geometry. 
The drift wave coupling effect is described by Vy acting upon yj y P — e0/T e and J\\, while 
interchange forcing is described by /C acting upon p e and [82|. In the case of tokamak 



edge turbulence, the drift wave effect is qualitatively more important 
important role for K, is to regulate the zonal flows 



80|, while the most 



831 ] . Detailed accounts on the role of 



magnetic field geometry shape in tokamaks and stellarators on plasma edge turbulence can 
be found in Refs. 84], |85| and [8Q, l87l|. in particular with respect to effects of magnetic field 
line shear and curvature js9| . 

An example for typical experimental parameters are those of the ASDEX Upgrade (AUG) 
edge pedestal plasmas in L mode near to the L-H transition for Deuterium ions with M, = 
Mr,: electron density n e = 3 • 10 13 cm -3 , temperatures T e = T{ = 70 eV, magnetic field 
strength B = 2.5 T, major radius R = 165 cm, perpendicular gradient length L± = 4.25 
cm, and safety factor q = 3.5. 

The dynamical character of the DALF/GEM system is determined by a set of param- 
eters characterising the relative role of dissipative, inertial and electromagnetic effects in 
addition to the driving by gradients of density and temperature. In particular, for the 
above experimental values, these are collisionality C = 0.51e(^ e Lj_/c s )(m e /Mj) = 5, mag- 
netic induction = e(47rp e / B 2 ) = 1, electron inertia fx = e(m e /Mi) = 5 and ion inertia 
6 = (qR/L ± ) 2 = 18350. 

The normalised values are similar in edge plasmas of other large tokamaks like JET. 
The parameters can be partia. 



University of Stuttgart 
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ly obtained even by smaller devices like the torsatron TJ-K at 
911 ]. which therefore provides ideal test situations for comparison 
between simulations and the experiment 92]. 

A review and introduction on drift wave theory in inhomogeneous magnetised plasmas 
has been presented by Horton in Ref. [93J], although its main emphasis is placed on linear 
dynamics. An excellent introduction and overview on turbulence in magnetised plasma and 



its nonlinear properties by Scott can be found in Ref. 94j, and a very detailed survey on 



drift wave theory with emphasis on the plasma edge is given by Scott in Refs. [95|, 196 ]. 



However, no tokamak edge turbulence simulation has yet reproduced the important 
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FIG. 6: Computation of plasma edge turbulence in the magnetic field of a divertor tokamak fusion 
experiment, using the DALF3 model as described in the text [Background figure: ASDEX Upgrade, 
Max-Planck-Institute for Plasma Physics]. 

threshold transition to the high confinement mode known from experimental fusion plasma 
operation. The possibility to obtain a confinement transition within first principle compu- 
tations of edge turbulence will have to be studied with models that at least include full 
temperature dynamics, realistic flux surface geometry, global profile evolution including the 
equilibrium, and a coupling of edge and SOL regions with realistic sheath boundary condi- 
tions. In addition the model still has to maintain sufficient grid resolution, grid deformation 
mitigation, and energy plus enstrophy conservation in the vortex/flow system. 

Such ,, integrated" fusion plasma turbulence simulation codes are currently under de- 
velopment. The necessary computing power to simulate the extended physics models and 
computation domains is going to be available within the next years. This may facilitate in- 
ternational activities (for example within the European Task Force on Integrated Tokamak 
Modelling) towards a ,, computational" tokamak plasma with a first-principles treatment of 
both transport and equilibrium across the whole cross section. The objective of this extensive 
project in Computational Plasma Physics is to provide the means for a direct comparison 
between our theoretical understanding with the emerging burning-plasma physics of the next 
large international fusion experiment ITER. 
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